Response time of a normal-superconductor hybrid system under the step-like pulse 
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The response of a quantum dot coupled with one normal lead and a superconductor lead driven 
by a step-like pulse bias Vl is studied using the non-equilibrium Green function method. In the 
linear pulse bias regime, the responses of the upwards and downwards bias are symmetric. In this 
regime the turn-on time and turn-off time are much slower than that of the normal system due to 
the Andreev reflection. On the other hand, for the large pulse bias Vl, the instantaneous current 
exhibits oscillatory behaviors with the frequency HQ = qVl- The turn on/off times are in (or shorter 
than) the scale of 1/Vt, so they are faster for the larger bias Vl- In addition, the responses for 
the upwards and downwards bias are asymmetric at large Vl- The turn-on time is larger than 
the turn-off time but the relaxation timei depends only on the coupling strength V and it is much 
smaller than the turn-on/off times for the large bias Vl- 

PACS numbers: 73.23.-b, 74.25.Fy, 74.50.+r 
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I. INTRODUCTION 

In the past two decades, nanoscopic physics has de- 
veloped significantly and becomes an active field of 
condensed-matter physics. The quantum transport prop- 
erty also becomes one of the most interesting phenomena 
in nanoscopic physics because of the possibility of de- 
signing and fabricating artificial setups in the nanometer 
scale. Based on the transport physics in nanoscopic sys- 
tem, a rich field for basic and applied research is opened* 2 - 
Furthermore, the time dependent nanoscopic transport, 
in which the external time dependent fields drives the 
electrons tunnel through a nanoscopic system, has re- 
ceived increasing attention in recent years. The main 
feature of the transport in the nanometer scale is that 
the electron keeps the phase coherence when traversing 
through the device. While the external time dependent 
field affects the phase factor of the incident electron dif- 
ferently in different parts of the system* 3 - If the external 
time dependent field is sinusoidal (e.g. microwave ra- 
diation), an electron can tunnel through the system by 
emitting or absorbing photons giving rise to the photon- 
assisted tunnelling (PAT). Electron transport with PAT 
has been extensively investigated for various systems, 
such as single or two coupled quantum dot (QD),— 
Kondo regime* 7 - hybrid system^ and so on. For transient 
transport, one of the most interesting issues is how fast 
can a device turn on or turn off a current. With the devel- 
opment of the molecular devices, there is clearly a need 
to technologically provide a particular viable switching 
device. Indeed, some recent experimental and theoreti- 
cal works have already begun to study the response of ac 
signals of the molecular devices.— Consequently, a step 
or pulsed ac signals are the simplest choice, since it can 
provide a less ambiguous measure of time scales. For 
this reason, the pulsed field was studied in a variety of 
systems, including Kondo regime ; 10 ' 11 a single QDji 2 . or 



nano structure , 13 ' 14 . 

So far, the study of response of pulsed bias is only 
focused on normal nanostructures. Since the interplay 
between nanoscopic physics and the physics of supercon- 
ductivity has made the hybrid structure a very fruitful re- 
search field,— it will be interesting to study the dynamic 
response a hybrid structure with a superconductor lead 
where the Andreev reflection is present near the normal- 
superconductor (N-S) interface. Indeed, there are many 
interesting phenomena in the N-S hybrid systems. First 
of all, because there exists an energy gap A in the su- 
perconductor, an incident electron from the normal side 
with energy e inside the gap A can not tunnel into the 
superconductor. But the tunneling can occur via a two- 
particle process, in which the incident electron is reflected 
as a hole with the energy — e. At the same time, a Cooper 
pair is created in the superconductor region. This is the 
Andreev reflection's Secondly, for the superconductor- 
normal region-superconductor (S-N-S) system, Andreev 
bound states form in normal region due to the Andreev 
reflections at N-S interfaces^ These bound states exist 
in pairs, and a Josephson supercurrent can flow through 
the S-N-S system which is carried by the Andreev bound 
stated Thirdly, when the S-N-S device is under an ex- 
ternal dc bias V , an ac current with frequency lu = 2\e\V 
appears. The time-average current versus bias V exhibit 
the subharmonic gap structure when eV < 2Aj^ 

In this paper, we explore the effect of Andreev reflec- 
tion on the ac response of hybrid system. Specifically, we 
investigate ac response of a quantum dot (QD) with a sin- 
gle level eo connected by a normal and a superconductor 
lead (N-QD-S). For simplicity, we consider a large QD 
so that the intra-dot electron-electron (e-e) is weak and 
can be neglected* 2 ^ The transient transport is driven by a 
pulsed bias potential W(t). For simplicity, the ac pulsed 
bias is only added in the left lead, and we set Wn(t) = 0. 
We consider two different pulsed bias: (i) upwards pulse 
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with W L (t) = for t < and W L {t) = V L otherwise, 
(ii) downwards pulse with Wl(€) — V L when t < and 
Wi(t) = otherwise. For normal structures, Wingreen 
et al. presented a general formula for the current driven 
by the time dependent external fields by using the non- 
equilibrium Green function (NEGF) method.^ With 
this general formula the time dependent current driven 
by the ac pulse can be calculated. For hybrid structures, 
the system is in steady state at t < and the current is 
time independent. At t = 0, bias is abruptly turned on 
for the upwards pulse case or turned off for the down- 
wards pulse case. After that, the system begins to relax 
and the Andreev reflection plays an important role in 
the relaxation process. Finally, the system enters into 
a new steady state. We find that, the relaxation time 
depends on the coupling strength and is slower in the N- 
QD-S system (named hybrid system hereafter) than in 
the N-QD-N system (named normal system hereafter). 
In the linear bias regime, the rising and falling processes 
are symmetric so that the turn-on time is same as the 
turn-off time. In this regime, the Andreev reflection is 
important. As a result, the instantaneous current shows 
a clear increase (decrease) before reaching the new steady 
state for the downwards (upwards) pulse. For the large 
bias case, the time dependent current oscillates with the 
frequency w = qVh- In this regime, the upwards and 
the downwards processes are asymmetric and the turn- 
on time is much larger than the turn-off time. In this 
nonlinear regime, the Andreev process is negligible and 
the current in the hybrid system is close to that of the 
normal system. 

The rest of this paper is organized as follows: In Sec. II, 
the theoretical formula for calculating the time depen- 
dent current in N-QD-S system is presented. To under- 
stand the numerical results, the current away from the 
current at t — is expanded to the first order in the 
external bias. In Sec. Ill, we show the numerical results 
along with some discussions. Finally, the brief summary 
is given in Sec. IV. 



II. THEORETICAL FORMULA 

Considering a hybrid system that consists of a QD cou- 
pled to a normal metal lead and a superconductor lead 
with the external time dependent bias potential Wt,{£) 
that is added only on the left normal lead. The Hamil- 
tonian of the system is written as follows: 



the following forms i 21 ' 22 



H = H T . + Hn + H t 



H 7 



(1) 



where Hl and H R describe the left normal lead and the 
right superconductor lead, respectively. Hp is Hamilto- 
nian of the isolated central QD, and the Ht couples the 
left and right leads to the QD. They can be written in 



H L 
Hr 

H d 
H t 



,k<r 



J2(eL,k + W L (t))Cl krT C L „ 

fe<7 

E £R,kC R k<J C RM + 
^r[AC Rj kiC R - k i + AC R _ k ^C Rkl ] 

k 

E ^A ha d. + h.c, (2) 



where a = L. R. The operator d a and C a ka destroy an 
electron with spin a in the QD and in the left or right 
lead, respectively. For simplicity, we only consider a sin- 
gle level in the QD and neglect intradot electron-electron 
Coulomb interaction. Under the adiabatic approxima- 
tion, the time-dependent bias potential can be included 
in the single electron energy eL,k(t). We separate ei,,k(t) 
into two parts: ei,,k and Wc{t), where th,k is the time- 
independent single electron energy and Wi,{t) is a time 
dependent part from the external time dependent bias 
potential. In this paper, WiS$) is the step- like pulse with 
two different forms: (i) upwards pulse with Wi,(t) = 
when t < and Wi(t) = Vl otherwise, (ii) downwards 
pulse with W L {t) = V L when t < and W L (t) = oth- 
erwise. These two types of pulse describe the system 
abruptly turned on or turned off at time t = 0. A in the 
Hamiltonian Hr is the superconducting energy gap. We 
assume that A is a real parameter by selecting a special 
phase of the superconductor lead in our calculation. 23 
Due to the existence of the superconducting lead, it is 
convenient to introduce the Nambu representation! 2 ^ In 
the Nambu representation, the Fermi energy of the left 
normal lead is set at the superconducting condensate and 
for the spin down electron the energy is negative and is 
viewed as the hole. So, the Hamiltonian in Eqs.© can 
be rewritten in the matrix form: 



t f e Ltk + W L (t) 
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(4) 



The current from the left lead to the QD can 
be calculated from the evolution of the number op- 
erator of the electrons in the left lead, N^^m = 
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T.k^ckm^aMd)^^ Usin S the Keldysh equation 
and the theorem of analytic continuation, the current 
through the left normal metal lead is expressed as: 3 ' 22 i 25 

J L {t) = -2qHe I dt'{[G r (t,t')Y<<{t',t) 

J — OO 

+G<(i,f )££(*', *)]n - [G r {t,t')Y><{t',t) 
+G<{t,t>)Y,l(t',t)] 22 } (5) 

Here the Green function G r l < and the self-energy 
are all two dimensional matrices in the Nambu represen- 
tation. Since the spin up and spin down are symmetric 
in the Hamiltonian, the current contributed by the elec- 
trons with spin up is same as the current by the spin 
down electrons. Consequently, the current is given by: 

J L (t) = -4gRe f dt'[G r (t,t')Y,<(t',t) 

J —OO 

+G<(t,t')El(t',t)} 11 . (6) 

Because of ££(*,*') = = {iT L /2)8{t-t')l (see 

Appendix) where I is the 2 unit matrix. Note that only 
G<(t,t) instead of G<(t,t') is needed in the Eq.©. By 
using the Keldysh equation G < — G r Y, < G a with the self- 
energies obtained in the appendix, the Green function 
G K {t,t) can be solved: 

G<(t,t) = fdh [ dt 2 G r (tM)^{hM)G a (t 2 ,t) 

a. J 

= i f ^/HG>)r fl ,HG>) + 
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f(w)A Li<r (tj, t)s a T L (u)Af(uj, t). 



(7) 



where a = ±1 denotes the spin up | and spin down |, 

r R ( |w| a 



T R (u) = 6{u - A) 



Vuj 2 - A 2 V A l w 



and 
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s l = 




1 



(8) 



where Bn = uo + e + iFl/2 + ii>T R [3/2, B 22 = uu — e + 
%T L /2 + ivT R p/2, p = A/Vu 2 - A 2 , 0' = oj/^/u 2 - A 2 , 
Det = B n B 22 + (r fl /3') 2 /4, and v = 1 for w > -A and 
1/ = — 1 otherwise. In the above derivation, the wide- 
band limit has been used and r Q are assumed indepen- 
dent of It also is worth mentioning that the Green 
function G r ^ a iuS) is not affected by the time-dependent 
bias potential Wl (t) ■ 

Substituting G < (t,t) [in Eq.([7|)] and the self-energies 
t) (in appendix) into Eq.©, the time-dependent 
current Jl(£) is obtained straightforwardly. Similar to 
the work in the normal system by Wingreen, Jauho, and 
Meir^ the current Jl(£) can also be split into two terms 
Jf(t) and J° L ut (t): 

J™(t) = AqJ %f{u)lm{T L [A L ^,t)]ii} 
J™\t) = -2qf^f(u)Re{r L {G r (u;)T R (u;)G a (u;) + 

ZaALaMs^LA+^U^jn}. (11) 

and J L (t) = J l L n {t) - Jl ut {t). Here the current J L n {t) is 
contributed by the electrons tunnelling from the left lead 
to the empty QD, and the current J£ u *(i) describes the 
electrons tunnelling from the QD to the empty left lead, 
so they have the opposite signi 3 - 

The above formulations [Eqs. (|9ll0lll|) ] for calculating 
the current are valid for any time-dependent bias Wl (t) ■ 
In the following, two special cases for upwards and down- 
wards pulses Wl{£) are substituted into these formula- 
tions to obtain AL<j{e,t) [Eq.@] and then the currents 
Jf{t) and J° L ut {t) [Eqs.d]}]. 

For the downwards pulse with < 0) = Vl and 

W L (t > 0) = 0, A L f(e,t) is found to be: 



A LZ5 , T (co,t<0) = G r (cu + V L ) 
A LDA (uj,t> 0) 



CT(u) + I ^e-^-^GUE) 



1 



1 



E-uj~V L -iO+ E — ui — i0 + 

(12) 



For the upward pulse with Wl (t < 0) = and Wl (t > 
0) = Vl, A L1 {e,t) is: 



A L<a (uj,t) 



dt 1 G r (t,t 1 )e Mt - tl)+iaI H dt2WL{t2 \(9) 



The Green functions G r ' a {uj) in Eq.jTJ) are the Fourier 
transformation of G r / a {t,t') with G r ' a {uj) = J d(t - 
t')e tuJ ( t - t ">G r/a (t,t'). Notice that in the present system 
the retarded and advanced Green functions G r ^ a (t, t') are 
still the function of the time difference t — t', although 
there exists the time dependent bias WL(t), since G r (ui) 
can be obtained from Dyson equation: 



Det 



B 11 ivT R p'/2 



ivT R p'/2 



Bo 



(10) 



A LUA (LO,t<0) = G r {uj) 
ALu, r {u,t >0) = G r {Lo + V L ) 

1 



dE 
2mi 



-i{E—u>—VL)tQr(jjr\ 



1 



E 



V L -i0+ E-lu- ifr 



(13) 



Here Ai CT (e, t) for the downward and upward pulse biases 
have been labelled by Ald,<j(€, i) and Alu,ct(£, t), respec- 
tively. For t < 0, the system is in the steady state, so 
Ald,<j(£, t) and A L u,a{ € : t) are independent of time t. On 
the other hand, for t > 0, they are obviously dependent 
on time t. For the purpose of numerical calculation, we 
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rewrite Ald/u.<j{^, t) for t > in the following form by 
using the residue theorem: 

A L D,i(u,t>0) = G r (u) + 

d T e i( - u+v ^ T G r (r) - / dre MT G r (r), 
A LUA (uj,t > 0) = & \lo + V l ) + 
e lVLt / dre lulT G r {T) - I dre l{u;+VL)T G r {t). 

(14) 

The expressions of AL,i(u>,t) are similar to that of 
ALf(u>,t) and can be obtained from Eq. (fl"4|) by chang- 
ing V L to — Vl- After solving G r {tu) and A La (u),t), the 
currents J L n (t) and J° L ut {t) [Eq.lQT])] can be calculated 
straightforwardly. In the limits t < and £ — > oo, the 
system is in the steady state. Ald.u^, t) and A £ [/ i0 -(o;, t) 
in Eq. (fT4j) then reduce to the value of the steady state in 
these two limits and so is the current Ji(t). For example, 
for the downward pulse, Aldm{^>, t) = G t (uj + <tVl) for 
t — ► , and Ald,<t{oj, t) = G r {u>) when t — ► oo. Further- 
more, the current Jz,(i) reduces to the one of the steady 
case with dc bias Vl when t < , and is zero when 
t — ► oo. On the other hand, for the upwards pulse, the 
current Ji(t) is zero when t < 0, and is same with the 
steady state current with the dc bias Vl in t — ► oo limit. 

In the small pulse bias Vl limits, we can expand 
ALo-(iL>,t > 0) to the first order of Vl as: AL a (w,t > 
0) = A £(y («,t - 0) + A\ a (u,t > 0)V L . A\ a {u,t > 0) 
can be expressed as: 

4>.'>«) = 
-iatj™ dre WT G r {T) - a f* drire^ G r \t) 

A 1 LU Ju,t>0) = 
iatj t °° dTe luJT G r (r) + a f* drire^ G r (r) (15) 

From Eq. (fT5"|) . we can see that A LDa (ui,t) — 
—A^ U(T (uj,t). This means that the upwards pulse and 
downwards pulse induce the same relaxation process in 
the small pulse bias Vl limits, except that the cur- 
rents deduced from them are relaxed in the opposite 
direction. Finally the currents J™(£) and J£"*(t) in 

small Vl limits can also be expanded as: J^ n *(t) = 
jrn/out^ + X iWo^(t)V L . Here X ln l out (t) is the first 
order expansion coefficient with the respect to Vl, and 
X m / out (t) is expressed as: 

X"\t)=4q J ^Tmf(u)r L {Ai tr {u,t)} n 
X<"*(t) = -2q J ^Re/Mr L ^ 

a 

{A\Au,t)s a Y L G a {u) + G\u) Sa YL[A\Au,t)\^ 

(16) 
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FIG. 1: (Color on line) The first order expansion coefficient 
X(t) of the current Jl(£) — Jl(0) vs. the time t for the down- 
wards and upwards pulse bias case in the hybrid N-QD-S sys- 
tem (a) and the normal N-QD-N system (b). The parameters 
are: r = 1, SF = 0, A = 15, e = 0. 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

In the numerical calculation, we set temperature to 
zero. In fact, finite temperature only makes the current 
curve more smooth and does not affect main features. 
We focus on the weak coupling case with Yl/ r <C A and 
set r = + Yr = 1 as energy unit. The energy gap of 
the superconductor is A = 15. The energy level eo in the 
central region is assumed to be zero which is same to the 
right Fermi level. Because at t < the system is in the 
steady state and the current is time independent, so we 
only plot the current Jz,(t) and the related quantities for 
t > in the following discussion. 

First of all, we study the small pulse bias Vl limit, in 
which the instantaneous current Jz,(i) can be expanded 
as: Ji(f) = Jl(0) + X(t)VL, and we also take the sym- 
metric barriers, i.e., 5Y = Yl — Yr = 0. The first- 
order expansion parameters x u ^ D ' m ^ out (t) of the cur- 
rents J L n (t) and versus the time t are plotted in 
Fig.l. Here the indices U and D denote the upwards 
and downwards pulses, respectively. For comparison, we 
also show the corresponding parameters x u I D > ln I out (t) 
for the normal system in Fig. 1(b). From Fig.l, we can 
see that the expanding parameters X(t) for the upwards 
and downwards pulses are symmetric, i.e. X u,m l out (t) = 
_ X D,in/out(ty It means t h at i n the small V L limit (i.e. 

the linear regime), the current turned off or turned on 
by the downwards or upwards pulses in exactly the same 
manner with the same time scale for both normal sys- 
tem and hybrid system. In other words, the case of the 
downwards pulse is the reversal process of the upwards 
pulse. So in the following, we use the upwards pulse as 
an example in the linear region. 

At time t < 0, the driving bias is zero for the upwards 
case. The system is in equilibrium state so the current 
JY is zero and J^' m and j^ out cancel to each other. At 
t = 0, the bias is abruptly switched on. At t > 0, the 
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bias Wl (t) is kept at Vl all along, the electrons with the 
energy in the bias window begin to traverse through QD. 
As the time t increases, J^' ln and j^- out deviate from 
the initial value (t = 0). A net current gradually in- 
creases and the device is gradually turned on. As a result, 
for the time t from to about 0.5(27r/r), X U ' in {t) and 
X U ' out (t) gradually increase (see Fig.l). This increasing 
process is almost the same for the normal system and 
the hybrid system. For the normal N-QD-N device, the 
relaxation process completes near the time t = 0.5(27r/r) 
and X U ' out (t) is the half of X U ' in (t) at large time. On 
the other hand, for the hybrid N-QD-S device, the be- 
havior of X ' m (t) is approximative^ the same as that of 
N-QD-N at large time, but X U ' out {t) begins to decrease 
when t > 0.5(27r/r), and it goes to zero at the end of 
the relaxation process. So the current J^(t — oo) for 
the N-QD-S device is twice as large as that of the N-QD- 
N device. We interpret these properties as follows. For 
the normal system, the fact that X m (t = oo) is twice of 
X out (t = oo) is because X m (t = oo) and X out (t = oo) 
are respectively contributed by the electrons tunnelling 
from the left lead into the empty QD and from the QD 
into the empty left lead with the electronic energy u> be- 
tween and Vl, and in this energy range the distribution 
of the left lead is /l(w) = 1 but the distribution in the 
QD is (f L (ui) + f R (u)))/2 = 1/2 for t = oo. While for 
the hybrid N-QD-S system, after the bias is turned on, 
the Andreev reflection begins to play a role. For J^' m , 
there is not much difference between the normal and hy- 
brid systems, since the electrons always tunnel from the 
left lead into the QD in both systems. But for J^' out , 
instead of reflecting electrons from QD into the left lead 
in normal system, the Andreev process reflects back the 
hole out of QD, which makes j^ out decrease. Note that 
Ta can be expressed as22: 

Ta ~ 64w 4 + (r 2 + ST 2 ) 2 ' ^ 

in the small bias limit (u> « 0) and 6F = 0, nearly all 
of the incoming electrons participate in the Andreev re- 
flection. Because of this, J^' out (t — oo) goes back to 
the initial (t=0) value. So X OMt (t) decreases to zero at 
t = oo. 

Next, we study the case of large pulse Vl- Fig. 2(c) 
and (d) depict the currents and J L n versus time t for 
the large pulse strength Vl = 10. For comparison, J£ ut 
and J l L n for the small pulse strength Vl — 0.1 are also 
plotted in the Fig. 2(a) and (b). The currents </£"* and 
J™ in the large bias case have the following characteris- 
tics: (i) In the small bias limit, the relaxation processes 
of upwards and downwards are symmetric. However, in 
the large pulsed bias Vl case, they are asymmetric (see 
Fig. 2c and 2d). For larger pulse bias Vl, the asymmetry 
are stronger, (ii) For the large bias case, J^' m for the 
upwards pulse oscillates with the frequency Ml = qVL, 
which can be clearly seen in Fig. 2c and 2d for Vl = 10. 
At Vl = 0.1 the oscillation disappears because Ml = qVL 
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FIG. 2: (Color on line) The currents J v J D > in and j^ D - out 
vs. the time t for the small pulse bias Vl =0.1 (upper panels 
(a) and (b)) and the large pulse bias Vl = 10 (lower panels 
(c) and (d)) in upwards and downwards pulse case. The left 
panels (a) and (c) are for the N-QD-N system and the right 
panels (b) and (d) are for the N-QD-S system. The other 
parameters are same with Fig.l 



is too small to oscillate before the system is completely re- 
laxed, (hi) J^' out (J®' out ) of hybrid system increases (de- 
creases) in the first and then decreases (increases), and it 
reaches maximum (minimum) before the current relaxed 
completely. This is different from the normal system, in 
which the currents J^ out and j®' out are monotonously 
relaxed into the steady state, (iv) The decreasing (in- 
creasing) process of the current j^- out (J^ ,ou *) in the 
large bias case is much weaker than that of the small bias 
case (see Fig. 2b and d). Because for the large pulse, the 
energy of the incident electrons w is large, then Ta -C 1 
from Eq. (fTT)) and the Andreev reflection is weak. So most 
of the incident electrons participate in the normal reflec- 
tion. Consequently, J£ K * is humped up (or down) slightly. 

Since the currents </£"* and J L n can not be observed 
independently, in the following we study the total current 
Ji(t) (Jl — J™ — J™ 1 ) which can be measured in the 
experiment. Fig. 3 shows the current J^' D driven by the 
upwards and downwards pulses versus the time t for the 
different pulse strengths Vl- Here the current responses 
to the upwards and the downwards pulse are symmetric 
at small linear bias Vl (see inset of Fig. 3), but are asym- 
metric at the large bias Vl (see main of Fig. 3). At the 
large Vl, Jl oscillates with the frequency Ml — Vl- On 
the other hand, j£ always changes slowly regardless of 
the large and small Vl- 

Now we focus the turn on/off time (or rise/fall timei^) 
and the relaxation time (or saturation time^i). The for- 
mer describes how fast can a device turn on/off a current, 
which is necessary to provide a particular viable switch- 
ing device, and the latter was referred to how fast can the 
device goes to a new steady state after a bias is abruptly 
switched on. For the small bias Vl, the turn-on time, 
turn-off time, and the relaxation time are almost same 
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Time (2it/r) 

FIG. 3: (Color on line) The currents Ji(t) vs. the time t 
for the N-QD-S system (a) and the N-QD-N system (b) with 
the different pulsed bias Vl- Main Figure is for the case of 
Vl = 1 and Vl ~ 10. The case of Vl = 0.1 is plotted in 
inset panel. The curves are labelled as: (1) Jl(Vl = 0.1); 
(2) J v l {Vl = 0.1); (3) j£(V L = 1.0); (4) J%(y L = 1.0); (5) 
Jl(Vl = 10); (6) Jl(Vl = 10). The other parameters are 
the same as Fig.l 



regardless of the normal and hybrid systems. However 
these (turn on/off or relaxation) times for the normal N- 
QD-N device are much shorter than that of the hybrid 
N-QD-S device. For the normal device, it has been well 
turned on or off at t — 0.2(27r/r). But for the hybrid de- 
vice, the system is turned on or off until t = 1.0(27r/r). 
On the other hand, for the large bias, the current J Lit) 
of the hybrid system has the same character with that 
of the normal system, so do the turn-on/off time and 
the relaxation time. Note that these three time scales 
are not equal now. The turn-on time is the fastest, even 
faster than the scale l/Vz,(27r/r). The turn-off time is 
in the scale 1/Vl(2tt/T), which is longer than the turn- 
on timeiiS The relaxation is ~ 0.5(27r/T), which is the 
longest and only depends on the coupling strength T. Let 
us explain why the character of J L (t) for the normal and 
hybrid system are the same at large Vl but very different 
at small Vl- Because at the large bias Vl, most of the 
incoming electron have the large energy lo, then Ta <C 1 
from Eq. (fT7|) and the Andreev reflection is weak, so the 
N-QD-S device and the N-QD-N device have the same 
turn-on/off and relaxation time. But for the small bias 
Vl, the resonant Andreev reflection is dominant in the 
transport process in the hybrid system, so that the cur- 
rent Ji(t = oo) of the hybrid system is twice as that of 
the normal system, and their character of Jl (i) also are 
very different. So we will only discuss the small pulsed 
bias Vl case further in the following. 

At last, we consider the case of asymmetric barriers 
(i.e., ST = Tl — Tfl ^ 0) and in the small pulsed bias 
Vl ■ Because in the small Vl the time-dependent current 
Jl (t) for the upwards and downwards pulse are symmet- 
ric, we only study the upwards case. Fig. 4 plots the cur- 



0.0 0.5 1.0 1.5 2.0 
Time (2?i/r) 

FIG. 4: (Color on line) The current JE(f) vs. the time t in 
the small pulse bias Vl = 0.1 with the different asymmetric 
coupling strength ST. The panel (a) and (b) are for the N-QD- 
S system and the N-QD-N system, respectively. The other 
parameters are same to the Fig.l 



rent J^it) versus the time t for the different asymmetric 
coupling strengths ST, and they have the following be- 
haviors: (i) As ST (i.e. Tl) increases, the current J^(t) 
rises faster, i.e. the turn-on time is shorter, because elec- 
trons with the energy in bias window can tunnel through 
the left barrier more easily with the larger T l ■ This rising 
process of Jj? are nearly same for the normal and hybrid 
systems, (ii) After the rise of (at t ~ 0.2(27r/r)), 
Andreev reflection begins to dominate and gives rise to 
different sequent relaxation processes for the normal and 
hybrid systems. At ST < 0, (t) of the hybrid system 
humps slightly in the relaxation process, which is obvi- 
ously different from the normal system in which J% (t) is 
monotonically relaxed into steady state. When ST = 0, 
J^(t) of the hybrid system passes a step and increases 
again. The relaxation time for the hybrid system is much 
longer than that of the normal system when ST < 0. 
When ST > 0, the relaxation processes of J^(i) are sim- 
ilar for the hybrid system and the normal system. These 
behaviors can be interpreted by combining the density of 
state (DOS) of the QD with the Andreev reflection pos- 
sibility T A . In fact, at ST > 0, the DOS of the QD in the 
hybrid system is similar to that of the normal system and 
Ta *C 1, so that the two systems have the similar turn- 
on/off and relaxation characteristic. On the other hand, 
when ST = or <5r < 0, the resonant or the near resonant 
Andreev reflection occurs, Andreev bound states appears 
in the QD, and the DOS of the QD is very different from 
the normal system. This makes the relaxation processes 
very different for the N-QD-S and N-QD-N systems, (iii) 
Although JY(t) for = +a and ST = —a (a is an 
arbitrary real number) experience different rising and re- 
laxation processes, they have the same steady value at 
t = oo. In fact, in the steady state case and at the small 
bias Vl limit, the transmission possibility of the normal 
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N-QD-N device is: 



T(u) 



r 2 - sr 2 



Auj 2 



p2 ' 



and the Andreev reflection possibility of the hybrid N- 
QD-S device is^ 



T A (u) 



(r 2 - «5r 2 ) 2 



4(4^ 2 + r<5r) 2 + (r 2 - sr 2 ) 2 ' 

with the current expressions Jl = —2q J — Vl) ~ 

f(u))T(u) and J L = -2qJ^(f(uj -V L ) - f(cu + 
Vl))Ta(w), respectively. Here T and Ta are the same 
for ±<5r when uj — 0, consequently J Lit = oo) also are 
same for ±&T. 



IV. CONCLUSIONS 

In summary, we have studied the dynamic response 
of current to the external upwards or downwards pulsed 
bias for the hybrid N-QD-S system. In the small bias 
Vl limit, the turn-on/off time and the relaxation process 
for the upwards and the downwards pulse bias are sym- 
metric. Comparing wtih the normal N-QD-N system, 
the Andreev reflection dominates the transport process. 
This makes the turn-on/off time much longer and new 
steady state current almost doubled. For the asymmetric 
barriers, the transport properties of the hybrid N-QD-S 
system are nearly same with the normal N-QD-N system 
when Tl > Tr. On the other hand, while Tl < Tr the 
current humps in the relaxation process which reflects 
the properties of the superconductor. Beyond the lin- 
ear bias regime, the rising process for upwards bias and 
the falling process for downwards bias become more and 
more asymmetric with the increasing bias Vl ■ The turn- 
on time is faster than the turn-off time, and the current 
versus the time t oscillates with the frequency Ml = Vl- 



ACKNOWLEDGMENTS 

This work was supported by the Chinese Academy of 
Sciences and NSF-China under Grant Nos. 90303016, 
10474125, and 10525418. J.W. is supported by RGC 
grant (HKU 7044/05P) from the government SAR of 
Hong Kong and LuXin Energy Group. 



APPENDIX 



In this appendix, we give the self-energy Y7' < for cou- 
pling to the left normal and right superconductor lead. 
To consider the wide-band limit, in which the hopping 
elements t^, a is independent with the momentum k and 
the density of state of the leads P^/ R (E) is energy inde- 
pendent, the self-energies from the coupling to the 
left normal lead with the time dependent bias potential 
Wz,(f) and in the Nambu representation are: 



k,L 

/Jk.Ldkt.L^'^^L 



k,L 



i I ^ / ( a ,)r ie --(t'-t)-/ t t '* 1 ^(tx) (19) 

y] t k,L9ki.,L( t ' '' ffifc,£ 

k,L 



i I ^(i_/( w ))r Le M*'-t)-H/* t '*iWL(*0 
J 2ir 

i f ^/( & ,)r £ e-M*'-*)-Hj?'<ft l w t (ti)(2o) 
J 2tt 



Here r £ = 27r|ifc £ | 2 p^ r , g^At' ,t) is the Green function 
of the isolated left lead, and f(ui) is the Fermi distri- 
bution. Notice that the retarded self-energy T, r La (t',t) 
is not affected by the time dependent bias Wi(t), so it 
is still the function of the time difference t' — t. Since 
the time dependent bias Wit) is applied only on the left 
normal lead and Wn(t) = 0, so the self-energies for cou- 
pling to the right superconductor lead are same with the 
steady state case and they can be written as&22 



2 V^ 2 



£<M = i0{w - A)/ B (w) 



uj A 
T R 1 



(21) 



2 V^ 2 - A 2 V A w 



uj A 



(22) 



where Tr = 27r|tfejj| 2 p^, A is the energy gap of the 
superconductor lead, and v = 1 for uj > —A and v = — 1 
otherwise. 
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